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ABSTRACT 

We analyze a 3D Smoothed Particle Hydrodynamics simulation of an evolving and later 
collapsing pre-stellar core. Using a 3D Continuum Radiative Transfer program, we generate 
images at 7, 15, 175 /im, and 1.3 mm for different evolutionary times and viewing angles. We 
discuss the observability of the properties of pre-stellar cores for the different wavelengths. For 
examples of non-symmetric fragments, it is shown that, misleadingly, the density profiles derived 
from a ID analysis of the corresponding images are consistent with ID core evolution models. We 
conclude that ID modeling based on column density interpretation of images does not produce 
reliable structural information and that multi-dimensional modeling is required. 

Subject headings: infrared radiation — ISM: clouds — ISM: dust, extinction — submillimeter — stars: 
formation 
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1. Introduction 

Molecular cloud cores are thought to be the di- 
rect progenitors of stars. However, their initial 
properties and early evolution are still poorly un- 
derstood. Current observations therefore aim to 
find and study these cores in more detail. Density 
structure, velocity field, and temperature distri- 
bution of the gas and dust are key parameters for 
the physical interpretation of the long lifetimes of 
these objects. The continuum radiation spectra 
of deeply embedded sources contain only ambigu- 
ous information about the density and tempera- 
ture distribution of the dust. Column densities 
can be inferred from the analysis of continuum im- 
ages in the mm-range (e.g. Ward- Thompson et al. 
1994) and the mid-infrared (MIR) (e.g. Bacmann 
et al. 2000). The emission of molecules within 
the core contains information about the struc- 
ture, velocity field (Tafalla et al. 2004, and ref- 
erences therein), and the turbulence (Ossenkopf, 
Klessen, & Heitsch 2001, and references therein). 
The currently discussed supporting mechanisms 
against gravitational collapse are magnetic fields 
and turbulence. There are a number of simula- 
tions treating the full 3D structure of collapsing 
pre-stellar cores (e.g. Bate, Bonnell, & Bromm 
2003; Ballesteros-Paredes, Klessen, & Vazquez- 
Semadeni 2003; Burkert & Bodcnheimer 1993, 
2000; Krumholz et al. 2003; Klessen, Heitsch, & 
Mac Low 2000; Heitsch, Mac Low, & Klessen 
2001). In some of these papers, the resulting col- 
umn density along a line of sight was compared 
with observationally obtained column densities. 
The question which structures of the distribution 
can be seen at which wavelength, however, can 
only be answered by producing images of the core 
for given dust properties. In turn, the majority 
of models that have been applied to observational 
data are mostly based on spherical symmetry (e.g. 
Andre, Ward-Thompson, & Motte 1996). They 
rely on spherically symmetric models of isolated 
star formation, which describe core formation, 
gravitational collapse, and proto-stellar accretion. 
Commonly, averaged radial density profiles are de- 
rived and compared with power-law density dis- 
tributions (e.g. Andre, Ward-Thompson, & Motte 
1996), or Bonnor-Ebert spheres (e.g. Alves 2004). 
The major difficulties using 3D models in observa- 
tions are i) the information loss due to the projec- 
tion effect, ii) the complex and unique structure of 



each individual pre-stellar cores, and iii) the nu- 
merical effort of multi-dimensional radiative trans- 
fer. 

In this letter, we investigate how an evolving 
cloud core simulated by a 3D Smoothed Particle 
Hydrodynamics (SPH) code would appear at dif- 
ferent wavelengths, which structures are visible, 
and which density distributions would be inferred 
using common ID models. The results from the 
SPH simulation are described in Sect. 2, along 
with a discussion of underlying assumptions. We 
present the images from the 3D Continuum Ra- 
diative Transfer (CRT) modeling and discuss the 
observability of different structures and physical 
effects. The comparison to density structures ob- 
tained from applying ID models to the images is 
given in Sect. 3, and the findings are summarized 
and discussed in Sect. 4. 

2. Cloud core evolution model and radia- 
tive transfer modeling 

2.1. 3D SPH simulation of the evolution 
of a cloud core 

We have calculated the evolution of a cloud core 
using a three-dimensional SPH code (version de- 
scribed in Bate, Bonnell, & Price 1995), originally 
developed by Benz, Cameron, Press, & Bowers 
(1990). The smoothing lengths of particles are 
variable in time and space, following the constraint 
that the number of neighbors for each particle has 
to be approximately constant with -/V nc i g h = 50. 
The SPH equations are integrated using a second- 
order Runge-Kutta-Fehlberg integrator with indi- 
vidual time steps for each particle (Bate, Bonnell, 
k Price 1995). 

The simulation was initiated with a mass of 
M = 3 Mq , adopting a spherically symmetric non- 
rotating homogeneous cloud with a temperature 
of T = 10 K, a diameter of d = 0.12 pc (corre- 
sponding to a density of 2 x 10~ 16 kg/m 3 ), and a 
mean molecular weight fx = 2.36 x 10~ 3 kg/mol. 
This configuration is Jeans unstable. A turbulent 
velocity field is added only at the beginning of 
the simulation with a Mach number of M. = 2 
and following an approximate Kolmogorov law 
P(k)dQ.k ~ fc~ 2 for the different modes. The tur- 
bulence supports the cloud core against collapse 
for the first 10 5 yrs. This enables the formation 
of a pre-stellar core- like structure, self-consistently 
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result of turbulent energy dissipation. A vari- 
able equation of state is used: isothermal for den- 
sities less than 10 16 molecules/m 3 and adiabatic 
for larger densities. 

Deviating from earlier work, the initial condi- 
tions are arranged in a way that the core reaches 
a dynamical equilibrium of density structures and 
velocity field before it evolves into a runaway col- 
lapse. The resulting structure is visualized by 
iso-density surfaces shown in the left panels of 
Fig. 1. The top left panel shows the early stage 
of core formation 5.6 x 10 4 yrs after the initial- 
ization (iso-density of 4 x 1CP 16 kg/m 3 ). Tur- 
bulence dominates the structure formation and 
creates several filamentary low-mass density max- 
ima. The duration of this period before the on- 
set of the collapse and thus the total "age" of 
the pre-stellar core stage depends on how turbu- 
lence is injected initially and its dissipation. In 
the course of time, the local density enhancements 
merge. After some additional 8.5 x 10 4 yrs just at 
the edge of gravitational instability, a single core 
has formed (second left panel, 5 x 10~ 17 kg/m 3 ). 
The kinetic pressure support breaks down due to 
rapid dissipation of turbulent energy inside the 
over-dense region and the core starts to collapse 
(t = 1.69 x 10 5 yrs, third left panel, 5 x lCT 17 
kg/m 3 ). A new single hydrostatic core forms when 
the gas becomes optically thick and the cooling 
time exceeds the dynamical time. Lateron, the 
central part of the core is replaced by a sink 
particle. In the bottom left panel, the struc- 
ture has flattened substantially, 20% percent of 
the total mass is already accreted onto the sink 
particle, and a massive disk has formed through 
an instability. It contains additional low mass 
condensations and independently, a second frag- 
ment has started to form with a hydrostatic core 
(5 x 10 -17 kg/m 3 ). The right panels give the iso- 
densities for 0.16, 0.5, 1.6, and 5.2xl0~ 18 kg/m 3 
for a time of 2.4 x 10 5 yrs, respectively. With in- 
creasing density, the second condensation becomes 
visible (see also animation at http://www.mpia- 
hd . mpg . dc /homes /stein/ Ani /animcf . htm) . 

2.2. 3D CRT modeling of the cores 

The SPH density distributions of the gas were 
discretized on a 3D grid and scaled to dust par- 
ticle distributions assuming a dust-gas mass ra- 
tio of 1/100 and an efficient gas-dust mixing. 



The dust number densities were processed with 
a 3D CRT code (Steinacker et al. 2003, 2002a; 
Stcinackcr, Bacmann, & Henning 2002b; Pas- 
cucci et al. 2004), producing 640 images of the 
cloud core at different wavelengths, times, and 
viewing angles, respectively. The temperatures 
were calculated from the radiative heating. Heat- 
ing by compression is irrelevant during the pre- 
stellar core phase due to the fact that the cool- 
ing timescale is much faster than the dynamical 
timescale. For the illustrative purpose of this let- 
ter, we used standard dust opacity data (Drain 
& Lee 1984) and a standard interstellar radia- 
tion field (Black 1994). Some of the images are 
shown in Fig. 2 for the wavelengths 7, 15, 175, and 
1300 /im (from top to bottom) and evolutionary 
times of 5.6 and 14.1 xlO 4 yrs, respectively (left 
to right). The wavelengths arc chosen to cover 
common observational windows (e.g. ISOCAM, 
ISOPHOT, IRAM, JCMT, CSO, SPITZER, HER- 
SCHEL). All images are scaled to have maximal 
contrast. A 10% random background noise repre- 
senting a mean background variation was added 
for illustrative purpose only. In the MIR, as ex- 
pected, the core is visible in absorption and the 
images reveal much of the outer thin structure es- 
pecially for the early stages. Detection of the in- 
ner, at later stages flattened structure is difficult 
and requires a careful background analysis. For 
wavelengths larger than 90 /mi, the cold dust can 
be seen in emission, revealing more of the inner 
structure at high densities, as the core also starts 
to get optically thin. This emission is dominating 
the mm images. 

Animations showing the images for all view- 
ing angles at different times of the evolution, 
as well as visualizations of the 3D density data 
cube can be found under http://www.mpia- 
hd. mpg.de/homes/stein/Ani/animcf.htm. 

3. ID analysis of the maps 

Projection effects are a severe source of mis- 
interpretation for structures seen in absorption 
or emission, as pointed out already, e.g., by 
Ballesteros-Paredes & Mac Low (2002). In Fig. 3, 
we show as an example two structures seen at 7 
/im. The upper left panel depicts an elongated fil- 
ament at an early stage of the evolution (5.6 x 10 
yrs), while the upper right panel gives a flattened 
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structure at a later stage (24.4 x 10 4 yrs). In the 
middle panels, the viewing angle was changed un- 
til we see the structures as a core-like feature. In 
the lower panel, they are zoomed and re-binned 
to an ISOCAM resolution typical for a core at 
150 pc distance. We determined the ID num- 
ber column density N(R) with the radius in the 
plane of the sky R by azimuthally averaging over 
annuli. As the absorption patterns have ellipti- 
cal shape, we have used elliptical annuli. The 
resulting number column density was inverted to 
a ID number density distribution n(r) with the 
radius r using recursive integration. In Fig. 4, 
we show the results for the early stage-filament in 
the main panel and for the later stage-disk in the 
inlet. The range of profiles n(r) which have been 
transformed from ./V(i?)-profiles along individual 
directions is given by the solid thin lines and the 
direction-averaged profile is plotted as solid thick 
line. The dashed line indicates the slope of a den- 
sity distribution following a r~ 2 -power-law as it 
was derived from ID core evolution models (Shu 
1977). Although these absorption maxima are 
slightly less extended than commonly observed 
cores, the ID model seems to provide a reasonable 
description of the derived distributions. It could 
be inferred from this fit that the underlying den- 
sity structure has an elliptical shape with a profile 
that - transformed to a spherical distribution - 
is in agreement with ID core evolution models. 
To compare with the true underlying 3D density 
distribution, we calculated the number density n 
for a grid of equally-sized cells from the SPH den- 
sity distribution. For each cell, we determined the 
distance to the center of the "core" -distribution 
defining a point in the n(r)-diagram. This point 
distribution was rcbinncd to a greyscale image 
for better clarity, where black refers to maximum 
number of cells per bin. The advantage of this 
representation is that a ID core with radial power 
law appears as a line in the n(r)-diagram with a 
gradient representing the powerlaw index. 

The true density distributions overlayed as 
grayscale-image are far from being lines as the 
filament is not ID spherically symmetric. The 
agreement with the ID evolution models would 
tend to validate static core formation models, al- 
though the core formation mechanism modeled 
here is highly dynamical. This is in agreement 
with the findings of Ballesteros-Paredes, Klessen, 



& Vazquez-Semadeni (2003). We also modelled 
the flattened structure at a later stage. The ra- 
dial profile within the disk will be visible in our 
grayscale-representation of n(r) as a line with the 
slope of the power-law exponent, and indeed the 
inlet in Fig. 4 shows a pronounced branch of the 
disk-like structure. If the core-flattening is con- 
firmed by an independent source of information, 
the column density profiles (and the resulting den- 
sity profiles) can be used to determine the radial 
profile of the disk-like structure. 

4. Conclusions 

We have presented 3D simulations assuming 
that initially low mass condensations pass through 
a stage of turbulence dominated condensation 
where they accumulate mass and merge together 
to form extended pre-stellar core like objects. 
The typical density structures in the cores are 
non-spherical throughout their evolution. The 
asymmetry is driven by the turbulent motion and 
causes complex structures from the very begin- 
ning. This complexity is partially seen in images 
that have been calculated from the densities ob- 
tained in the cloud core simulation. However, 
projection effects can lead to a severe misinterpre- 
tation of images. We showed that a ID analysis 
of the vicinity of the density maxima would sug- 
gest a density profiles in agreement with ID core 
collapse models. The underlying density struc- 
ture, however, is intrinsically three-dimensional 
and deviates strongly from the obtained ID model 
distribution. 

As the column density also enters the opti- 
cal thin emission in the mm-range (aside from 
the Planck function), we expect the same pro- 
jection ambiguities to occur when interpretating 
mm-maps of dense molecular cloud regions. This 
aspect will be discussed in a forth-coming paper. 

We conclude that ID modeling based on col- 
umn density interpretation of images does not pro- 
duce reliable structural information. For flattened 
structures appearing in later stages of the core evo- 
lution, a 2D modeling might be applicable, but 
for the general case, multi-dimensional continuum 
and line radiative transfer modeling is required to 
derive consistent density and temperature distri- 
butions of the gas and dust in pre-stellar cores. 
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Fig. 1. — Iso-density surfaces of the averaged SPH 
density distributions of a cloud core fragment. 
The left panel shows densities of 40, 5, 5, and 5 
xlCT 17 kg/m 3 at the times 5.6, 14.1, 16.9, and 
27.2 xlO 4 yrs after the start of simulation, respec- 
tively. The right panel gives iso-density surfaces 
at the time 24.4xl0 4 yrs and densities of 0.16, 0.5, 
1.6, and 5.2xl0 -18 kg/m 3 , respectively. 
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Fig. 2. — Images of the cloud core fragment at 
the wavelengths 7, 15, 175, and 1300 /xm (left 
columns from top to bottom), and at the times 
5.6 and 14.1 x 10 4 yrs after start of simulation (left 
to right), respectively. 




Fig. 3. — Examples of features in 7 (ira images 
mimicking dense cores. The top panels show struc- 
tures which look like a core when choosing an 
appropriate viewing angle (middle panels). The 
lower panels zoom into the core-like structures and 
switch to ISOCAM resolution. 
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Fig. 4. ID density profiles n(r) obtained from 
the 7 /xm-image by azimuthally integrating the 
column density around the absorption maximum 
along elliptical tracks. The thin solid lines mark 
the range of profiles for individual directions, and 
the direction-averaged profile is represented by the 
thick solid line. The dashed line indicates an r~ 2 - 
dependency. The true 3D density distribution 
usph was discretized on a cell grid, and trans- 
formed to a point distribution in the n(r)-plane. 
The number of points is shown as gray-scale image 
where black means maximum number of grids with 
a certain density. The main picture corresponds 
to an elongated filament from the early stage, and 
the inlet to a flattened structure from a later stage 
of the core evolution, respectively. 
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